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We present a theoretical analysis of the effects of uniaxial magnetic anisotropy and contact-induced 
exchange field on the underscreened Kondo effect in S = 1 magnetic quantum dots coupled to 
ferromagnetic leads. First, by using the second-order perturbation theory we show that the coupling 
to spin-polarized electrode results in an effective exchange field -B c ff and an effective magnetic 
anisotropy D e ff- Second, we confirm these findings by using the numerical renormalization group 
method, which is employed to study the dependence of the quantum dot spectral functions, as well 
as quantum dot spin, on various parameters of the system. We show that the underscreened Kondo 
effect is generally suppressed due to the presence of effective exchange field and can be restored 
by tuning the anisotropy constant, when \D c h\ = |-B e ff|- The Kondo effect can also be restored 
by sweeping an external magnetic field, and the restoration occurs twice in a single sweep. From 
the distance between the restored Kondo resonances one can extract the information about both 
the exchange field and the effective anisotropy. Finally, we calculate the temperature dependence of 
linear conductance for the parameters where the Kondo effect is restored and show that the restored 
Kondo resonances display a universal scaling of S = 1/2 Kondo effect. 

PACS numbers: 75.75.-c,75.76.+j,72.15.Qm,85.75.-d 



I. INTRODUCTION 

Although manifestation of the Kondo effect in 
nanoscopic systems of spin S > 1/2 has been the sub- 
ject of extensive experimental and theoretical studies for 
more than a decadep^ it is still attracting considerable 
attention. From the experimental point of view, this was 
triggered by a rapid development of techniques^ allow- 
ing for controlled preparation and investigation of single 
magnetic impurities, such as atoms and molecules, placed 
on a surface^TJI] or captured in a junctionFSHin] Further- 
more, an important issue is the interaction of individ- 
ual large-spin atoms or molecules with the environment, 
which may contribute to a magnetic anisotropyf^HED A 
significant uniaxial magnetic anisotropy, in turn, results 
in an energy barrier for switching the molecules's spin be- 
tween two metastable states - the feature indispensable 
for pot ential applications in information storage tech- 
nologies IMlll Interestingly enough, the magnetic state of 
such a system can in pri nciple be controlled by means 
of spin-polarized cmrentsP™2 which has already been 
experimentally confirmed™ 

In order to be able to exploit advantageous features 
stemming from the presence of magnetic anisotropy, a 
possibility of its external control would be very desirable. 
Indeed, several experiments have so far confirmed the 
feasibility of such a control. The most straightforward 
way to modify the magnetic anisotropy of an adatom 
is just to change its nearest atomic environment, which 
can be achieved simply by deposition of t he adato m at 
topologically different points of a substrate! 25 ! 34 ! 35 ! More 
elaborate techniques demonstr ated for molecules involve 
the application of electric field 1121221 or even the mechan- 



ical modification of the molecular symmetryPS In fact, 
the latter method allows for a fully-controllable and con- 
tinuous tuning of the anisotropy constant, which was 
demonstrated for a spin S = 1 quantum dot in the un- 
derscreened Kondo regime. EHHU 

Screening of a quantum dot spin appears when the 
dot becomes strongly coupled to electrodes. For temper- 
atures T smaller than the Kondo temperature Tk, the 
spin exchange processes due to electronic correlations can 
lead to an additional sharp peak in the density of states 
- the Kondo- Abrikosov-Suhl resonance. Generally, in 
order to observe full screening of a magnetic impurity 
spin S, t he imp urity should be coupled to 2S screening 
channels! 37 ! 44 ! 45 ! In turn, a typical experimental setup for 
measuring transport through quantum dots or molecules 
involves usually two contacts. This implies that when 
connecting the spin S = 1 dot to two (say first and sec- 
ond) leads, the spin could be in principle fully screened.^ 
In order to observe the underscreened Kondo effect, one 
needs to use a more specific setup, as demonstrated by 
Roch et alW^ Since the screening becomes effective when 
T < and the Kondo temperature depends exponen- 
tially on the dot-lead coupling strength T, by connecting 
the dot asymmetrically to external leads one obtains two 

1 ( 2) 

different Kondo temperatures: T K for the first (second) 
lead. The underscreened Kondo effect can be then ob- 
served when the condition ^ « T « T| is fulfilled^ 
In such a case, the spin is only partially screened by elec- 
trons of the strongly-coupled lead, while the other lead 
serves as a weakly coupled probe. Despite its theoret- 
ical simplicity, the first experimental realization of the 
under screen ed Kondo effect was reported only very re- 
centlyPEH 
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In this paper, motivated e.g. by the experiments 
of Parks et oZ.J^ we analyze the transport properties 
of a spin 5 = 1 system strongly coupled to a fer- 
romagnetic reservoir. In particular, we focus on dis- 
cussing how the uniaxial magnetic anisotropy and the 
ferromagnetic-contact-induced exchange field affect the 
underscreened Kondo effect. Our analysis is based on 
the full density-matrix numerical renormalization group 
(fDM-NRG) methodj^lEIl which is known as the most 
powerful and exact in addressing transport properties 
of various nanostructures in the Kondo regime. As we 
are mainly interested in the aspects of the underscreened 
Kondo effect, which are related to the coexistence of 
magnetic anisotropy and ferromagnetism of the screening 
channel, we assume a model in which only one electrode 
is attached to the dot, see Fig. [T] Such a set up defines 
a typical one-channel Kondo experiment j 36 l 46 l where the 
role of the second (weakly coupled) electrode in the for- 
mation of the Kondo resonance can be neglected (the 
corresponding Kondo temperature tends to zero). Nev- 
ertheless, the second electrode, being a weakly coupled 
probe (e.g. a tip of STM), will enable the measurements 
of the conductance through the system and the local den- 
sity of states. 

The paper is organized as follows. In Sec. II we de- 
scribe the model and method used in calculations. Sec- 
tion III is devoted to basic concepts, where we discuss 
the spectrum of an isolated dot and the effects of renor- 
malization due to the coupling to electrode. Numerical 
results and their discussion are presented in Sec. IV with 
a special focus on the effects due to magnetic anisotropy. 
Finally, the conclusions can be found in Sec. V. 

II. THEORETICAL DESCRIPTION 
A. Model 

The total Hamiltonian of a two-level magnetic quan- 
tum dot coupled to an external lead, see Fig. [T] can be 
written as 

7~L = T~L mo \ + Hlcad + 7~Lt\in, (1) 

where the first term describes the quantum dot, the sec- 
ond one refers to the lead, whereas the final term repre- 
sents tunneling processes between the dot and the lead. 

A bare two-level quantum dot can be characterized by 
the model Hamiltonian 

■Hmol = E E ^ n 3 a + U '^2 n 3t n jl + U ' E n ^ n ^' 
j=l,2tr=t,4. 2=1,2 o-,o-'=t4 

+ JS 2 + DS 2 Z +B Z S Z . (2) 

In the above, n,j a — Cj (3 Cj a and c^c^) denotes the cre- 
ation (annihilation) operator of an electron with spin a 
and energy Sj in the jth level (j — 1,2). For conve- 
nience, we write the energy levels as ei = e — 5/2 and 




two level quantum dot 

FIG. 1. (Color online) Schematic of a two-level magnetic 
quantum dot coupled to a metallic ferromagnetic electrode. 
The magnetic moment of the electrode (denoted by a bold 
arrow) is collinear with the dot's easy axis. The dot levels 
have energies £i and £2, respectively, with S being the level 
spacing, while J denotes the exchange interaction. The tunnel 
matrix elements between the dot and the lead are denoted by 
Tj a for the dot level j and spin a. 

£2 = £ + (5/2, where £ is the average value of the two 
levels while 8 is the level spacing. The Coulomb en- 
ergy of two electrons of opposite spins occupying the 
same level is given by U (assumed to be the same for 
both levels) , whereas the inter-level Coulomb correlations 
are described by U'. For simplicity, we assume equal 
correlation energies U = U' in the following. Further- 
more, J stands for the interlevel exchange interaction 
with S = S\ + S2, where Sj is the electron spin oper- 
ator for the dot level j, Sj = § e ia <J aa >c ja , , with 
(7 = (<r x ,a v ,a z ) being the Pauli spin operator. Accord- 
ing to the Hund's rules, this interaction should be gener- 
ally of & ferromagnetic type (J < 0), nonetheless the pos- 
sibility of a weakly antiferromagnetic coupling has also 
been reported.^! Finally, the lowest order uniaxial mag- 
netic anisotropy is represented by the anisotropy constant 
D, and the last term of Eq. ([2| describes the Zeeman en- 
ergy of the dot in an external magnetic field B z applied 
along the dot's easy axis, with gfiB = 1. 

The Hamiltonian for a ferromagnetic metallic reservoir 
of noninteracting itinerant electrons is given by 

Hload = E E £ kv a la a ka- (3) 

k a=t,l 

Here, a} k(T {a kry ) creates (annihilates) an electron of energy 
£fc CT , where k indicates a wave vector, while a is a spin 
index of the electron. It is important to mention that 
in the following discussion we assume that magnetic mo- 
ment of the electrode remains collinear to the dot's easy 
axis. 

Finally, tunneling of electrons between the electrode 
and the dot is described in general by 

W tun = E E E T jkAa^ + H.C., (4) 
k j=l,2 CT =t4 

where Tj^ denotes the tunnel matrix element between 
the dot's jth level and the electrode, see Fig. [l] In the 
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following we assume that both levels are coupled sym- 
metrically to the electrode, i.e. T\k a = T^ka = 2fc CT -Al- 
though such foundation is not the most general one] 53 * 54 ! 
it is sufficient for the present analysis of the effects result- 
ing from magnetic anisotropy and exchange field in the 
context of the underscreened Kondo problem. In order 
to further facilitate calculations, we assume that the full 
spin-depend ence is included exclusively via the matrix 
elements T a where the fc-dependence has also been 
neglected. In addition, we assume a symmetric and flat 
conduction band extending within the range [— W, W], so 
that the density of states is p(oj) = p = 1/(2W), and we 
use W = 1 as the energy unit. Consequently, the spin- 
dependent hybridization function reads, T a = 7rp\T a \ 2 . 
Now, introducing the spin polarization coefficient P of 
the electrode, defined as P = (T t - r i )/(r T + TJ, 
the spin-dependent coupling can be parameterized as: 

r m) = r(i ± p), with r = (r t + rj/2. 

B. Objectives and method of calculations 

The main quantity we are interested in is the zero- 
temperature, spin-dependent equilibrium spectral func- 
tion of the quantum dot 



s 7 =o 



1 



(j,j' = l,2), (5) 



with (( c j a \ c l'a))ij standing for the Fourier trans- 
form of the retarded Green's function ((c J - cr |ct, tr ))j = 
— i8(t)({cj a (t), Cj, a (0)}). It is worth noting that the spec- 
tral function with two identical indices, i.e. A 3 J(uj), is 
related to the spin-resolved density of states associated 
with the jth level, whereas A>J (to) with j ^ j' cor- 
responds to processes of electrons entering and leaving 
the dot at different levels. Because measuring the spin- 
resolved components of the spectral function may pose 
a serious experimental challenge, we will focus on dis- 
cussing the total spectral function. Thus, we introduce 
the normalized full spectral function A(lo), 



A( u ) = f^^r^'( u ). 



(6) 
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The importance of the spectral function A{uS) stems 
from the fact that in a two-terminal setup with the 
second lead being a weakly-coupled probe, e.g. a tip 
of an STM microscope, the differential conductance of 
the system at bias voltage eV can be approximated as 
|^ ~ ^A(uj = el/)PM10n the other hand, the spectral 
function for uj — > 0, ^4(0), determines the linear-response 
conductance. 

In the light of the preceding discussion, the cen- 
tral problem is the calculation of the spectral function 
A J J (ui). In the Kondo regime, this can be reliably 
done b y mean s of the numerical renormalization group 
(NRG)p 4 * 48 * 49 ! which enables us to analyze the static 
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FIG. 2. (Color online) Sketch showing how the uniaxial 
magnetic anisotropy lifts partially the degeneracy between the 
components of the triplet S = 1. Note that in reality the 
energy of the S z = component is independent of magnetic 
anisotropy D. 



and dynamic properties of the system in the most ac- 
curate manner. The essential idea of the method lies in 
a logarithmic discretization of the conduction band and 
mapping of the system's Hamiltonian onto a semi-infinite 
chain, with the quantum dot residing at the initial site. 
Iterative diagonalization of the Hamiltonian by adding 
consecutive sites of the chain allows then for resolving 
key properties of the system at energy scale A - ™/ 2 , with 
A > 1 denoting the discretization parameter and n a 
given iteration. 

In order to address the present problem efficiently, the 
calculations have been performed with the use of the 
flexible density- matrix numerical renormalization group 
(DM-NRG) code.^Eal In this study we exploited the 
^charge (1) x f^spin(l) symmetries corresponding to con- 
servation of the electron number (charge) and the zth 
component of the total spin. 



III. BASIC CONCEPTS 
A. Isolated quantum dot 

Before presenting and discussing numerical results on 
the spectral functions, it is advisable to have a closer 
look at the energy spectrum of an isolated quantum dot. 
For the sake of clarity of the following discussion, let us 
assume that there is no external magnetic field, B z = 0, 
so that the system's behavior is entirely determined by 
both the sign and magnitude of the uniaxial anisotropy 
constant D, as shown schematically in Fig. [2j 

In order to observe the underscreened Kondo effect, 
the quantum dot needs to be occupied by two elec- 
trons which are ferromagnetically exchange-coupled. The 
ground state of the dot is then a triplet with the compo- 
nents, 



|T + ) = |5, = +l) = |t)i r t>2, 
\T ) = \S, = 0) = ~ [ 

\T_) = \S z = -l) = \i) 1 \i) 2 , 



|t)i|4) 2 + |4)i|t) ; 



(7) 



where \x)j denotes the local state of the jth level, with 
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X = 0, 4-j tj d corresponding to zero, spin-down, spin- 
up and two electrons occupying the level, respectively. 
As long as an external magnetic field and the magnetic 
anisotropy are absent, the three triplet states remain de- 
generate and st + — £t = £ t_ = 2e + U + 2 J, see Fig. [2] 
Moreover, the triplet remains the ground state provided 
the condition, 5/2 - 2U + 5J/4 < e < -5/2 - U - 5J/4, 
is satisfied.^ However, the magnetic anisotropy D lifts 
this degeneracy and the triplet becomes partially split, 



(a) 



(b) 



(c) 



e T+ = e Tl = 2e + U + 2 J 4 D, 

e To =2e + U + 2 J, 

e T _ = e Tl = 2e + U + 2 J 4 D. 



(8) 



As one can see, the energy of states \T + ) and |T_) de- 
pends on D, while the energy of |Tq) is independent of D. 
Consequently, when D < 0, the ground state is two-fold 
degenerate and corresponds to the states \T + ) and |T_), 
while for D > 0, the ground state corresponds to |T ), see 
Fig. [2} This will have a large impact on the Kondo effect, 
as discussed later on. Note that the presence of magnetic 
field additionally splits the states \T + ) and |T_). 



B. Effective exchange field and anisotropy 

The next question that straightforwardly arises is what 
happens when the S = 1 quantum dot becomes attached 
to a reservoir of electrons. If the temperature is lower 
than the Kondo temperature, T < Tk, the conduction 
electrons can then screen only a half of the dot's spin, 
whereas the residual spin-one-half is left unscreened. At 
zero temperature, the system behaves then as a singu- 
lar Fermi liqu id, i.e . the Fermi liquid with a decoupled 
S = | object.^^ In addition, it turns out that mul- 
tiple spin-flip processes responsible for the Kondo reso- 
nance lead to renormalization of the quantum dot pa- 
rameters. One can in principle distinguish two different 
effects associated with such a renormalization. First, the 
spin degeneracy of the dot i s lifted by an effective tunnel- 
induced exchange field -B e ff- 55 1 56 1 ^ "^ It was shown that 
the exchange field can be tuned by a gate voltage^ 
and can be comp ensated by applying an external mag- 
netic fieldP^HSa The second effect, on the other hand, 
is related to the renormalization of the anisotropy con- 
stant, AD, which results in an effective anisotropy D e g, 
D c ff = D + AD, as shown schematically in Figs.[3ja)-(c). 
This renormalization is independent of the spin polariza- 
tion P of the lead, while it depends weakly on the gate 
voltage and cannot be compensated by external magnetic 
field. 

The renormalized energies (i = 0, ±) of the triplet 
state can be found from the second-order perturbation 
theory in the tunneling Hamiltonian as eV ; = + Se^ 
(i = 0, ±), where is the second-order correction 

of the respective triplet component. The renormalized 
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FIG. 3. (Color online) (a)-(c) Schematic representing level 
renormalization due to the effective anisotropy constant D c g 
and effective exchange field B e g for the quantum dot with 
5=1. Dashed lines in (a) and (c) correspond to energy 
levels without the renormalization of the uniaxial magnetic 
anisotropy constant, D e s = D, see Eq. (111. It is assumed 
that B e ff < 0, so that the ground state (in (b) and (c)) is 
|T+). Bottom panel: the dependence of B e ff and D e g on the 
anisotropy constant D for e/U = —1.4 (d) and on the level 
position £ in the Coulomb blockade regime for |-D|/T = 0.1 
(e). Note that in (e) for D < one also gets D e fi < 0, but for 
practical reasons \D e s\ is plotted. For comparison, in (d) we 
also show AD and D. The other parameters are: e = — 17. 5r, 
5 = 2.5r, U = 12. 5r, J = -2.5r, and P = 0.5. 



triplet energies may be written in the following way 

£t + = £t q + 5et q + AjS + B c ff + B z , 
e To = e To + 5e To , (9) 
e T _ = eTg 4. 5s Tg 4. £) cff _ B cS — B z , 

where the effective exchange field B c r induced by a fer- 
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romagnetic contact is given b; 



B, 



off : 




i - /M 



3=± 



(10) 



and the effective anisotropy D e g can be expressed as 

D off = D + AD (11) 

with the renormalization of the anisotropy constant AD 
of the form 



AD 



1 - f(u) 1 - f(u) 




(12) 



The prime superscript in the above equations symbol- 
izes Cauchy's principal value integrals, and f(w) stands 
for the Fermi-Dirac distribution function of the con- 
tact. We note that terms involving 1 — /(w) represent 
here electron-like charge fluctuations, due to which the 
molecule loses one electron, while terms with f(ui) re- 
fer to the hole-like processes, when the charge of the 
molecule is increased by one electron. Furthermore, 
E a ,p — £a ~ £p is the energy difference between the 
corresponding states. The respective energies of singly 
occupied states are, Si± = e ± 5/2 + 3 J/4 + D/4, while 
the energies of states with three electrons are given by, 
e 3± = 3e + 3(7 ± 5/2 + 3J/4 + D/4. In Eq. ([o]) 5e To 
(SeTo < 0) denotes the second-order energy correction of 
the triplet component |T ) (uniform shift of the whole 
triplet). The explicit form of 5et is not relevant for 
the present discussion, since it is the difference between 
the above energies of triplet components that determines 
the occurrence and features of the Kondo effect. In the 

IVC]is in the 



low temperature regime, which except Sec. 



main scope of the work, Eqs. (10) and (12 1 simplify sig- 
nificantly, 
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Et ,1- Etq,1+ Et ,3- Et ,3A 



(13) 



(14) 



We also note that the energies of triplet, Eq. ([9|, explic- 
itly include the external magnetic field B z . This will be 
relevant for the discussion of system transport properties 
in the presence of magnetic field, that will be presented 
in the next section. 



From Eqs. ( 10 ) and ([131 follows that the exchange field 



is an intrinsic effect resulting from the spin-dependence 
of tunneling processes, and vanishes for P — > 0. More- 
over, B c ff displays monotonic dependence on the level 
position e: with B c g — at the particle- hole symmetry 
point of the model, i.e. for e = —3U/2, and B e g ^ for 



e ^ — 3U/2. This is shown in Fig.[3]je). In addition, B c h 
also depends on the anisotropy constant D, see Fig.J^d). 
This dependence, however, is much weaker than the de- 
pendence on e. 

On the other hand, since ej\ = Et + D, one can imme- 
diately conclude from Eqs. ( 12 ) and (14 1 that AD — > 
as D — > 0. It is interesting to note that tunneling of elec- 
trons leads to suppression of the magnetic anisotropy, 
AD ^ for D ^ 0, in the whole Coulomb blockade 
regime. Furthermore, unlike the effective exchange field, 
AD is an even function of the level position e with the cx- 
tremum (maximum for D > and minimum for D < 0) 
in the particle- hole symmetry point, where 



AD, 



7T 



(2C/-5J-3D) 5 



4<S 2 



(2U - 5 J + D) 2 — A5 2 



(15) 



The dependence of the effective anisotropy D e g on the 
magnetic anisotropy constant as well as on the level po- 
sition is shown in Fig. |3jd,e). One can note that D c g 
depends strongly on D and only weakly on the level 
position, which is just opposite to the behavior of the 
exchange field B c g. Moreover, while B g is due to the 
spin-dependence of tunneling processes and vanishes for 
nonmagnetic leads, D c fj does not depend on the spin po- 
larization and is finite also when P = - as long as 
D ^ 0. 

The above discussion suggests that transport proper- 
ties should be mainly determined by the interplay of the 
effective anisotropy D c s, contact-induced exchange field 
.B e ff and the Kondo temperature Tr-. Additionally, the 
behavior of the total spectral function also significantly 
depends on the tunnel-coupling strength T (to observe 
the Kondo physics the coupling should be sufficiently 
large). Experimental results show that the Kondo phe- 
nomena in quantum dots can be observed when T is of the 
order of a few tent hs of meV for temperatures of the or- 
der of mK. I21512U22I Accordingly, in numerical calculations 
we assume r = 0.5 meV and use T as the relevant energy 
scale. For the quantum dot we assume the parameters 
that are comparable to those observed in experiments, 
e = -17.5r (e = -8.75 meV), 5 = 2.5T (5 = 1.25 meV), 
U = 12. 5r (U = 6.25 meV), and J = -2.ST (J = -1.25 
meV), if not stated otherwise. Note that we assumed 
e/U > —3U/2, so that if the lead is ferromagnetic, the 
ground state is \T + ) due to B e g < 0, see also Fig.^c). 



IV. NUMERICAL RESULTS AND DISCUSSION 

In the following we present and discuss numerical re- 
sults on the dot's spectral density as a function of the 
anisotropy constant D, spin polarization of the lead 
P, and external magnetic field B z . The main focus, 
however, will be on the effects arising from the mag- 
netic anisotropy and effective exchange field. Gener- 
ally, the magnetic anisotropy in systems under consider- 
ation can take fairly large values, and can range approxi- 
mately from |D| < 0.05 meV for single-molecule magnets, 
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| ii | i4 | i7 | 7 3 ] U p fa a f ew me y f or magnetic adatoms like Mn, 
Fe, Co, l 23 l 25 l 35 l or some magnetic moleculesP3 Further- 
more, Park et al— have shown that mechanical strain in 
a certain type of Co complexes allows for a fully control- 
lable and continuous change of the magnetic anisotropy 
of a molecule. This effect occurs since the stretching or 
squeezing of a molecule leads to modification of the crys- 
tal field exerted on the central Co ion. For the above 
reasons, the following results will be presented for a wide 
range of both positive and negative uniaxial anisotropy 
constant D. 



A. Influence of uniaxial magnetic anisotropy 

1. Nonmagnetic lead 

To begin with, let us first discuss briefly how the spec- 
tral function of the system depends on the magnetic 
anisotropy in the case of a nonmagnetic electrode!^ Fig- 
ure Ha) shows the dependence of A(0) on D for several 
values of the exchange interaction J. It can be seen that 
A(0) « 2 for \D\ < T K , where T K is the Kondo tem- 
perature, with A(0) decreasing below its unitary value 
once \D\ > Tk- Moreover, the resonance dies away 
more abruptly for D > 0, where the spectral function 
is practically equal to zero above some threshold value 
of the anisotropy constant. Accordingly, one should ex- 
pect there a vanishingly small linear conductance of the 
system. Indeed, such a behavior has been observed by 
Parks et aLp^who reported splitting of the Kondo peak 
due to stretching the molecule. The asymmetry between 
the decrease of ^4(0) for positive and negative D is asso- 
ciated with different ground states of the quantum dot, 
see Fig. 2] For D > and D c g > T K , the ground state is 
\Tq) and no spin-flip processes are possible, consequently 
A(0) becomes abruptly suppressed. On the other hand, 
for D < and \D c g\ > Tk, the ground state is two-fold 
degenerate, with equally occupied states |T + ) and |T_), 
see Eq. Due to the spin selection rules for tunnel- 
ing processes, the second-order spin-flip cotunneling is 
then suppressed and the Kondo resonance becomes sup- 
pressed as well [A(0) starts decreasing]. However, there 
are fourth-order tunneling processes that are still possible 
and therefore A(0) decreases rather slowly with increas- 
ing \D\, opposite to the case of positive D. 

The features discussed above depend on the Kondo 
temperature Tr-. Since Tr- is a function of the energy 
difference between the ground state and single and three- 
particle virtual states, the Kondo temperature can be 
tuned by changing the exchange interaction J. When 
increasing |J|, one effectively increases the energy dif- 
ferences and Tk becomes decreased. Accordingly, the 
suppression of the Kondo effect occurs for smaller values 
of D. This can be clearly seen in Fig. ga). Note that 
the width of the maximum in A(0) as a function of D is 
roughly equal to 2Tk- 
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FIG. 4. (Color online) Total normalized spectral function 
A(uj) of the dot for uj — > 0, ^4(0), as a function of the uniax- 
ial magnetic anisotropy constant D for several values of the 
interlevel exchange interaction J in the presence of (a) non- 
magnetic (NM) and (b) ferromagnetic (FM) electrode. The 
insets in (a) and (b) display ^4(0) as a function of D/Y plotted 
in a linear scale in the case of nonmagnetic (a) and ferromag- 
netic (b) electrode for different values of J. Note that in the 
inset in (a) the Kondo peaks have been shifted by —0.01 for 
better visibility, (c) Variation of the spectral function with 
the energy ui and uniaxial magnetic anisotropy constant D 
for J/r = —2.5. (d) Cross-sections of the plot in (c) for se- 
lected values of the uniaxial anisotropy constant D > 0. The 
other parameters are: e = — 17. 5r, S — 2.5F, U = 12. 5r, and 
P = 0.5 (except (a) where P = 0). 



2. Ferromagnetic lead 

The situation becomes much more interesting when the 
nonmagnetic reservoir is replaced by the ferromagnetic 
one. The dependence of the spectral function A(Q) on 
D is presented in Fig. W[b). First, for small values of 
\D\, and thus also |-D e ffTi the height of the Kondo res- 
onance is significantly reduced as compared to the case 
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FIG. 5. (Color online) (a) Total normalized spectral func- 
tion A(uj) of the dot for uj — > and (b) the corresponding 
expectation value of the dot's spin zih component as func- 
tions of the uniaxial magnetic anisotropy constant D. Except 
J/r = —2.5, remaining parameters are the same as in Fig. [4] 



of a nonmagnetic electrode, which is due to the presence 
of exchange field i? c ff. Second, as D > increases, one 
observes the revival of the Kondo effect at some resonant 
value of the magnetic anisotropy constant, D = D rcs . 
However, further increase in D results in the drop of 
the spectral function to zero, so that the behavior of the 
system for large and positive D resembles that of a dot 
coupled to a nonmagnetic electrode. In addition, the de- 
pendence on the exchange coupling is also qualitatively 
similar to that in the nonmagnetic case. When increasing 
| J\, Tk is reduced and the width of the Kondo resonance 
as a function of D becomes decreased as well, see the 
inset in Fig. |4^b). In fact, the width of the Kondo reso- 
nance for a given value of J is of the same order in the 
case of nonmagnetic and ferromagnetic leads, compare 
the insets in Figs. |4|a,b). In addition, the value of D res 
decreases with increasing J, which is due to the corre- 
sponding dependence of B c g and D c g on the exchange 
coupling J. 

Since the spectral function can be substantially modi- 
fied upon altering the anisotropy constant D, the follow- 
ing discussion will be focused on the interplay of D c g and 
-B ff, that governs the transport behavior in the Kondo 
regime. In the remaining part of the paper we will present 
and discuss numerical results for a fixed value of the ex- 
change coupling, J/r = —2.5 [corresponding to the bold 
lines in Figs. |4|a,b)]. 

As follows from Eq. (13 1, the strength of B c g depends 
on the spin polarization P of the reservoir, B c g ~ PT. 
Moreover, through the energy differences between re- 
spective states, B c g is a function of all parameters of 



the model, including the magnetic anisotropy constant 
D. First of all, unlike AD [see Eq. (Ill])], B c g is finite 
for D — > 0. More specifically, for parameters used in 
Figs. [3jd) and |4|c,d) it approaches a constant value of 
Beff/r « —0.087. Since the exchange field lowers the 
energy of the highest-weight triplet component \T+), sec 
schema (b) in Fig. [3j one can observe almost full spin po- 
larization of the dot, (S z ) — > 1, see Fig. [gb) for D -> 0. 
Note, however, that the dot's spin can be flipped to 
(S z ) -> -1 for e < -3U/2, where B cS > and |T_) 
becomes the ground state of the system. This can be 
achieved for instance by applying a gate voltage. In ad- 
dition, B e g also depends on the anisotropy constant and 
it can either increase or decrease depending on the sign of 
D, see the solid line in Fig.|3jd). For the parameters used 
in calculations, the modification of B c s for |D| = T/10 
is however rather small (~ 1%). The variation of the 
exchange field as a function of D is thus of rather minor 
significance for T ^ \D\, but nevertheless the interplay 
of D (D c g) and B e g turns out to be crucial for the oc- 
currence of the Kondo effect. 

Let us focus first on the case of D < 0, where no 
restoration of the Kondo resonance takes place, see the 
left side of Fig. |4jb) . The ground state for P = would 
be doubly degenerate. Because of the exchange field, 
this degeneracy, however, is lifted and the ground state 
is |T + ). The state |T + ) remains the ground state in the 
whole range of D < considered in this paper. In con- 
sequence, there is no Kondo effect for D < 0. 

The situation, however, is much more complex in the 
case of D > 0, where the restoration of the Kondo reso- 
nance appears, see the right side of Fig. |4^b). For small 
values of D, where the dominant energy scale due to 
renormalization processes is set by the exchange field 
-B c ff, the ground state is |T+) and the situation is sim- 
ilar to that for D < 0. Nonetheless, as the magnetic 
anisotropy grows, the condition B c g + D c g = becomes 
satisfied at some point (note that for the assumed pa- 
rameters B e g < while D c s > 0) and the state |T + ) 
gets degenerate with the state |To). The difference be- 
tween the spin zth components of the states \T + ) and 
[To) is 1, so the second-order spin-flip cotunneling pro- 
cesses are possible and the Kondo effect can be restored. 
A(0) reaches then its maximal value. 

The above described behavior can be also observed in 
the full energy dependence of the spectral function A(uj), 
see Figs. gc,d). The Kondo resonance is restored when 
D = D ICS ss r/10 and is immediately suppressed once 
D > D TCS . On the other hand, the suppression is less ef- 
fective on the left side of the restored Kondo resonance, as 
already discussed above. In addition, the energy depen- 
dence of the spectral function in Fig. (4^d) reveals small 
side peaks for D > D rcsi which are reminiscent of the 
Kondo effect and occur for energies oj ~ ±|-D ff + -B e ff| 
corresponding to restored degeneracy of the states |T + ) 
and | To)- Apart from this, at large energies, lj ~ ±C7, 
there are typical Hubbard resonance peaks. 

Since the occurrence of the Kondo resonance depends 
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Thus, the estimations based on the analytical expressions 
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FIG. 6. (Color online) (a) Dependence of the magnetic 
anisotropy constant D Yes at which the Kondo resonance is 
restored on the interlevel exchange coupling J for P = 0.5 
and (b) on the spin polarization P for J/r = —2.5. Other 
parameters as in Fig. [4] 



on the ratio of B c g and D e g, it is interesting to study 
variation of A(Q) with D for different values of lead's 
spin polarization. This is shown in Fig. |5ja). When 
P = 0, A(0) shows a maximum for such D that the con- 
dition |-Deff| ^5 Tk is fulfilled. If the spin polarization 
is finite, the maximum is shifted towards larger values 
of anisotropy and occurs precisely when the states |T + ) 
and \Tq) become degenerate, i.e. for D = D ICS . Note, 
that the width of the peak in ^4(0) as a function of D 
is of the same order for all values of spin polarization 
P. Since the dot is coupled to one electron reservoir, 
only half of the dot's spin can be screened by the con- 
duction electrons for D = D ICS . As a result, in the un- 
derscreened Kondo regime the expectation value of the 
dot's spin should reach (S z ) ~ 1/2, since the unscreened 
residual spin- 1/2 is polarized due to the presence of ex- 
change field. This can be seen in Fig. [5^b), which shows 
the dependence of (S z ) on D for several values of P. For 
D < D ICS , the ground state of the dot is \T + ) and (S z ) is 
close to one, while for D > D les , the ground state is |T ) 
and (S z ) = 0. On the other hand, for D rj D rcs , one finds 
(S z ) — 1/2, see Fig. [5jb) . However, closer analysis of 
(S z ) shows that it actually fails in attaining its maximum 
value for D < D TCS . This is related with the fact that the 
ratio U/T is relatively large for the assumed parameters 
and there is nonzero occupation probability of other spin 
components of the triplet.^ However, when increasing 
the spin polarization P of the lead, the splitting of the 
levels grows due to the exchange field, B c g ~ PT, and 
the occupation of the triplet component \T + ) is raised. 
In consequence, one finds that (S z ) — > 1, if P — > 1. 

Knowing the analytical condition, B e g + D e g = 0, for 
the occurrence of the Kondo resonance for D > 0, it 
is instructive to analyze the role of magnetic anisotropy 
renormalization. For this purpose, in Fig. [6] we present 
the dependence of D rcB on J and P. The solid (dashed) 
line corresponds to D res determined from the analytical 
formulas for B c g and D c g with (without) including AD, 
while the dots show D rcB obtained from NRG data. As 
one can see, the NRG results are in very good agreement 



(14), are quite satisfactory. The renormalization of D is 



thus an important effect that needs to be included in the- 
oretical considerations of spin S > 1 systems exhibiting 
magnetic anisotropy. 

To demonstrate additional features of the interplay be- 
tween magnetic anisotropy and exchange field, we show 
in Fig. [7] the energy and spin polarization dependence of 
the normalized spectral function A(uS). The full energy 
dependence of the spectral function may prove to be use- 
ful in predicting some qualitative information concerning 
transport properties of the system at a finite bias. As it 
was discussed earlier, the spin polarization P determines 
the strength of i? c ft, without affecting D c ff- In conse- 
quence, all w-dependent features in Fig. [7] should in prin- 
ciple stem from the changes of B e g with respect to D e g. 
For P < P rcs , where P rcs is the value of spin polarization 
at which the Kondo resonance is restored, one observes 
a well-pronounced dip, which indicates that the system's 
ground state is nonmagnetic, i.e. \T ). In the present pic- 
ture, increasing u) turns to be equivalent (to some extent) 
to the application of an external bias voltage when the 
dot would be asymmetrically attached to two contacts.^ 
For P — > 0, the spectral function attains then its local 
maximum at u> that approximately corresponds to the de- 
generacy of the states |To) and \T + ). From this, in turn, 
D e ff can be straightforwardly obtained, i.e. D c g u>. As 
P increases, the magnitude of the exchange field \B g\ in- 
creases as well, and this is accompanied by a decrease in 
the energy gap B c g + D c g. This appears then as a grad- 
ual narrowing of the dip, until the two states, i.e. |T + ) 
and | To), become degenerate, and the Kondo resonance 
is restored, see Figs. [7Fa)-(e). For larger values of P, the 
system's ground state is \T + ) and the Kondo resonance is 
suppressed again. There are however two satellite peaks 
at energies uj sa ±|_D e g + £? e ff|; whose position depends 
linearly on P, and whose height diminishes as P grows 
further. Moreover, it turns out that for larger values of 
D, see Figs. [7](f )- (g) , no restoration of the Kondo effect 
is possible. This is because the magnitude of exchange 
field is too low to compensate the effective anisotropy 
D e ff and the degeneracy of states cannot be restored even 
if P4 1. 



B. Influence of an external magnetic field 

Figure [8] shows the dependence of the spectral function 
A(0) and the average value of dot's spin (S z ) on the mag- 
netic field B z . It can be noticed that while the Kondo 
effect for D < can be restored just for a single value of 
B Zl for D > the restoration occurs twice. In order to 
understand this behavior one should bear in mind that 
the external magnetic field affects the components \T + ) 
and \T_) of the triplet state and thus can be used to 
compensate the splitting induced by exchange field B c g 
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FIG. 7. (Color online) Total normalized spectral function A(uu) of the dot shown as a function of energy uj and the spin 
polarization parameter P for indicated values of the uniaxial magnetic anisotropy constant D. Other parameters are the same 
as in Fig. [5] 



due to ferromagnetic contact,^ see Eq. (JoJ) . In the case 
of vanishing D and B z = 0, the Kondo resonance is sup- 
pressed due to the exchange field B e g and the system 
ground state is \T + ). With increasing B z , all the three 
components of the triplet state become degenerate once 
£? c ff + B z ?s 0, and the system is in the underscreened 
Kondo regime. Nevertheless, because only a half of the 
dot's spin can be screened by conduction electrons, the 
remaining spin-1/2 can be polarized by any infinitesimal 
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FIG. 8. (Color online) Dependence of the total normalized 
spectral function A(uj) for uj — > (a,b) and the average zth 
component of the spin (S z ) (c,d) on an external magnetic field 
B z oriented along the dot's easy axis. Different lines corre- 
spond to selected values of the uniaxial magnetic anisotropy 
constant D, and the left panel corresponds to D < 0, while 
the right panel refers to D > 0. Except for P = 0.5 and 
J/r = —2.5, all remaining parameters are the same as in 
Fig-H 



magnetic field (at zero temperature) . This leads to strong 
sensitivity of the ground state on magnetic field, which 
hinders the full restoration of the underscreened Kondo 
effect P 

For D < 0, the increase of B z can only restore the 
degeneracy between the highest and lowest-weight com- 
ponents of the triplet. In consequence, one observes a 
small peak at B z sa \B c g\, see Fig.^a), where the ground 
state changes from |T + ) to |T!_), see Fig. (8^c) . Note that 
the position of this peak does not depend on D. On 
the other hand, for positive anisotropy D > 0, the full 
restoration of the Kondo effect is possible, see Fig.[8^b). 
Moreover, contrary to single-level quantum dotsp^the 
restoration with increasing B z occurs twice. This can 
be understood by studying the evolution of the ground 
state with magnetic field, see Fig. Isfd). For B z = and 
D > 0, the ground state is a singlet, |T ). By lower- 
ing the magnetic field, B z < 0, the ground state changes 
to |T + ) for B z ps — D ff — B c s, while by increasing B z , 
once B z ss D e g — B c ff, the ground state changes to |T_). 
Consequently, once B z ps ±£W — B c g, the two-fold de- 
generacy of the ground state becomes restored and the 
Kondo effect can develop. One observes then two max- 
ima in A(0), see Fig. [8jb). Note, however, that different 
states are responsible for these two Kondo peaks. For 
B z ps — D c g — B e ff, it is the degeneracy between the states 
\T + ) and |To) that results in the formation of the Kondo 
effect, while for B z ss D cS — B cS , the states |T_) and 
\Tq) are degenerate. It is also worth noting that since 
the Kondo temperature is rather independent of D, the 
width of the restored Kondo peaks is the same for all 
values of D, see Fig.[8][b). 

From the magnetic field dependence of ^4(0) one can 
obtain the information about the magnitude of both B c g 
and D c g. Suppose the restoration of the Kondo ef- 
fect occurs for B z = B^l and B z = -Brci, then the 
effective anisotropy constant can be related to a half 
of the distance between the two restored Kondo reso- 
nances D eS = \b£ ] s - b£1\/2- On the other hand, 
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FIG. 9. (Color online) The temperature dependence of the 
normalized linear conductance G(T)/G(0) (a) and the uni- 
versal scaling curves (b). The relevant curves correspond to 

(i) underscreened Kondo effect, with D = 0, P = 0, B z = 0, 

(ii) restored Kondo resonance when B c g + D c h — 0, with 
D/V « 0.1, P = 0.5, B z = 0, and (iii) restored Kondo reso- 
nance by magnetic field B z = D c a — B ff, with D/V ~ 0.2, 
P — 0.5, B z /V w 0.22. The curves (ii) and (iii) display scaling 
typical for spin S =1/2 Kondo effect [dotted-dashed line in 
(b)] . Other parameters the same as in Fig.gwith J/T = -2.5. 



the magnitude of the exchange field can be found from 
the average, B cS = + B&})/2. Studying the 

magnetic field dependence of the zero-energy spectral 
function, which would correspond to measuring the low- 
temperature zero-bias conductance, may be thus useful in 
obtaining information about both the effective anisotropy 
and exchange field. 



C. Temperature dependence of linear conductance 

It is very instructive to study the temperature T de- 
pendence of the linear conductance, G{T), in the un- 
derscreened Kondo regime and for parameters where the 
restoration of the Kondo effect appears. Generally, the 
conductance in the underscreened Kondo regime can be 
measured by attaching a second weakly- coupled elec- 
trode, which due to much smaller K ondo t emperature, is 
irrelevant for screening the dot's spin BSEal The linear con- 
ductance has been calculated by means of NRG method 
with the full density-matrix, and the Meir-Wingreen for- 
mulaP The normalized linear conductance as a func- 
tion of temperature T is shown in Fig. [9][a) for the un- 
derscreened Kondo effect (dashed line), i.e. for P = 



D = B z = 0, and for parameters where the restora- 
tion of the Kondo effect occurs - first when the con- 
dition D c ff + B g = is met (solid line) and second 
when the restoration is obtained by applying magnetic 
field, i.e. when D c s — B c s — B z = is satisfied (dotted 
line) . The relevant Kondo temperatures are also given in 
the figure. Tr is defined here as the value of T where 
G(T)/G(0) = 1/2. Figure ^b) displays the universal 
scaling curves of normalized conductance G{T) /G(0) as 
a function of T/Tk- For D = 0, P — and in the absence 
of magnetic field, we observe scaling typical for the under- 
screened Kondo regime.^ which has been recently mea- 
sured experimentally.^ The temperature dependence of 
the linear conductance for parameters where the restora- 
tion of the Kondo effect occurs also turns out to be uni- 
versal, see Fig. 9^b), however the scaling is completely 
different from the underscreened Kondo effect. For pa- 
rameters where the restoration occurs, the ground state 
is two-fold degenerate, i.e. either \T + ) and |T ) or |T ) 
and |T_) components of the triplet state are degenerate, 
therefore one should expect the same scaling as for the 
spin 5=1/2 Kondo effect. ^Indeed, we compare the 
universal scaling of the linear conductance for restored 
Kondo resonances of S = 1 quantum dot with the scal- 
ing for typical S = 1/2 Kondo effect and find perfect 
agreement, see Fig. |^b). 



V. SUMMARY AND CONCLUSIONS 

By means of numerical renormalization group method, 
we have studied transport properties of a magnetic S = 1 
quantum dot coupled to a ferromagnetic lead in the un- 
derscreened Kondo regime. Due to the coupling of the 
dot to an external lead, the following two effective pa- 
rameters are shown to play important role: the effective 
exchange field B e g and the effective anisotropy constant 
-D ff. The interplay of the corresponding interactions is 
crucial to understand behavior of the system transport 
properties, especially regarding the evolution (suppres- 
sion or restoration) of the Kondo effect as a function of 
various parameters of the model considered. Using the 
second-order perturbation theory, we have derived an- 
alytical formulas for both B e g and D c g. It turns out 
that the effective anisotropy D e g depends strongly on 
D and D c ff — > as D — > 0. Furthermore, D c s is an 
even and weakly changing function of the level position 
e, with an extremum at the particle-hole symmetry point, 
e = — 3C//2, and does not depend on the spin polarization 
of the ferromagnetic lead. The effective exchange field 
-B c ff, on the other hand, depends linearly on lead's spin 
polarization P and B c s — > for P — > 0. Furthermore, 
it is an odd function of level position e and vanishes at 
the particle-hole symmetry point, e = — 3U/2. B c s also 
depends on D, although this dependence is rather weak. 
We compared the analytical formulas for B e s and D e g 
with the NRG data and found very good agreement. 

By performing extensive NRG calculations, we have 
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studied variation of the spectral functions with various 
parameters of the system. We have shown that the un- 
derscreened Kondo effect is generally suppressed due to 
the presence of magnetic anisotropy and exchange field. 
It can be, however, restored by tuning the magnetic 
anisotropy constant D. The restoration occurs only for 
positive anisotropy, D > 0, while no restoration takes 
place when the magnetic anisotropy is negative, D < 0. 
Moreover, the restoration of the Kondo resonance also 
occurs as a function of magnetic field applied along the 
easy axis. By sweeping the magnetic field, the Kondo ef- 
fect can be restored twice in a single sweep. The restora- 
tion always occurs due to the degeneracy between the 
components of the triplet state that differ in the spin S z 
quantum number by 1. 

We have also determined the temperature dependence 
of the linear conductance for some characteristic param- 
eters, where the restoration of the Kondo effect occurs. 
It turned out that the restored Kondo resonances exhibit 
a universal scaling as a function of T/Tk characteristic 
of spin 5 = 1/2 Kondo quantum dots. This is due to 
the fact that for parameters where the restoration of the 
Kondo effect is possible, the ground state is two-fold de- 
generate. 

Finally, we would like to emphasize that when con- 
sidering spin-resolved transport through nanostructures 



of spin S > 1 exhibiting magnetic anisotropy, there are 
two relevant and distinct effects that need to be taken 
into account in order to fully understand behavior of the 
system. The first one is the exchange field induced by 
ferromagnetic contact, and the second one is associated 
with effective (renormalized) magnetic anisotropy. We 
also remark that nanoscopic systems for which the mag- 
netic anisotropy is a generic feature, as the ones discussed 
in this paper, present just one possible way of employing 
magnetic anisotropy as a key element of novel spintron- 
ics devices. More recently, it has been suggested that 
spin-anisotropy can also be generated in spin-i sotrop ic 
systems by spin-dependent transport of electrons! 76 ! 77 ! 
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